HEAVY 
FLAVOURS^ 



Heavy Flavours 8, Southampton, UK, 1999 



PROCEEDINGS 



Determination of the weak phase 7 = ai'gCK*/,) 



Matthias Neubert 

Newman Laboratory of Nuclear Studies 
Cornell University, Ithaca, NY 14853, U SA 
E-mail: 



neubert Omail . Ins . Cornell . edi: 






o 

(N 



Abstract: Various strategies for extracting or constraining the weak phase 7 with controlled the- 
oretical uncertainties are reviewed. Measurements of the rates for the hadronic decays B^ -^ ttK 
provide largely model-independent information on 7. Hadronic uncertainties enter only at the level 
of nonfactorizable contributions to the decay amplitudes that are power-suppressed in A/mi, and, si- 
multaneously, violate SU(3) flavor symmetry or are doubly Cabibbo suppressed. In addition, these 
decays have a rich potential for probing physics beyond the Standard Model. Alternative proposals 
for learning 7 are also discussed briefly. 
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THE MAIN OBJECTIVES of the B factories are 
to explore the physics of CP violation, to de- 
termine the flavor parameters of the electroweak 
theory, and to probe for signs of physics beyond 
the Standard Model. This will test the Cabibbo- 
Kobayashi-Maskawa (CKM) mechanism, which 
predicts that all CP violation results from a sin- 
gle complex phase in the quark mixing matrix. 
Facing the announcement of evidence for a CP 
asymmetry in the decays B -^ J/ijj Kg by the 
CDF Collaboration [Q , the confirmation of direct 
CP violation vaK -^ tttt decays by the KTeV and 
NA48 groups [g[ || , and the successful start of the 
asymmetric B factories at SLAC and KEK, the 
year 1999 has been an important step in achiev- 
ing this goal. 

The determination of the sides and angles of 
the "unitarity triangle" V*^Vud+V*^Vcd+VlVtd = 
plays a central role in the _B-factory program. 
Adopting the standard phase conventions for the 
CKM matrix, only the two smallest elements in 
this relation, V*^^ and Vtd, have nonvanishing ima- 
ginary parts (to an excellent approximation). In 
the Standard Model the angle /? = —arg{Vtd) can 
be determined in a theoretically clean way by 
measuring the mixing-induced CP asymmetry in 
the decays B — > J/iJ^Ks- The preliminary CDF 
result implies sin 2/3 = 0-79toil 0- The angle 
7 = arg(y^*j), or equivalently the combination 



a = 180° — (3 — "f, is much harder to determine 
B. Recently, there has been significant progress 
in the theoretical understanding of the hadronic 
decays B — > ttK, and methods have been de- 
veloped to extract information on 7 from rate 
measurements for these processes. Here we dis- 
cuss the charged modes B — > ttK, which from 
a theoretical perspective are particularly clean. 
We will, however, also survey alternative meth- 
ods for learning 7 from other decays. 

In the Standard Model the main contribu- 
tions to the decay amplitudes for the rare pro- 
cesses B -^ ttK are due to the penguin-induced 
flavor-changing neutral current (FCNC) transi- 
tions b — *■ sqq, which exceed a small, Cabibbo- 
suppressed b — > uus contribution from W^-boson 
exchange. The weak phase 7 enters through the 
interference of these two ( "penguin" and "tree" ) 
contributions. Because of a fortunate interplay of 
isospin, Fierz and flavor symmetries, the theoret- 
ical description of the charged modes B^ — > nK 
is very clean despite the fact that these are ex- 
clusive nonleptonic decays [^Iq, Q . Without any 
dynamical assumption, the hadronic uncertain- 
ties in the description of the interference terms 
relevant to the determination of 7 are of rel- 
ative magnitude O(A^) or 0(esu(3)/^c), where 
A — sin 0(7 ~ 0.22 is a measure of Cabibbo sup- 
pression, esu(3) "^ 20% is the typical size of SU(3) 
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breaking, and the factor 1/Nc indicates that the 
corresponding terms vanish in the factorization 
approximation. Factorizable SU(3) breaking can 
be accounted for in a straightforward way. 

Recently, the accuracy of this description has 
been further improved when it was shown that 
nonleptonic B decays into two hght mesons, such 
as -B — > ttK and B — > tttt, admit a heavy-quark 
expansion g. To leading order in A/rrib, but 
to all orders in perturbation theory, the decay 
amplitudes for these processes can be calculated 
from first principles without recourse to pheno- 
menological models. The QCD factorization the- 
orem proved in [pj improves upon the phenome- 
nological approach of "generalized factorization" 
g, which emerges as the leading term in the 
heavy-quark limit. With the help of this theorem 
the irreducible theoretical uncertainties in the de- 
scription of the B —> ttK decay amplitudes can 
be reduced by an extra factor of 0(A/mf,), ren- 
dering their analysis essentially model indepen- 
dent. As a consequence of this fact, and because 
they are dominated by FCNC transitions, the de- 
cays B^ — > ttK offer a sensitive probe to physics 
beyond the Standard Model @, |l|, 0, |I|, |l|, 
much in the same way as the "classical" FCNC 
processes B -^ Xgj or B ^ Xs l'^l~ . 

In the following three sections we discuss how, 
in the context of the Standard Model, B^ — > ttK 
rate measurements can be used to derive bounds 
on 7 and to extract its value with controlled theo- 
retical uncertainties. We then explore how these 
analyses could be affected by New Physics. Fi- 
nally, we briefly summarize alternative methods 
to determine 7 using other decay modes. 

1. Theory of B^ — ^ nK decays 

The hadronic decays B — > 'kK are mediated by 
a low-energy effective weak Hamiltonian, whose 
operators allow for three distinct classes of fla- 
vor topologies: QCD penguins, trees, and elec- 
troweak penguins. In the Standard Model the 
weak couplings associated with these topologies 
are known. From the measured branching ratios 
one can deduce that QCD penguins dominate the 
B —> ttK decay amplitudes |Q, whereas trees 
and electroweak penguins are subleading and of 
a similar strength pa. The theoretical descrip- 



tion of the two charged modes B^ 



--K° and 



B^ — > TT^K^ exploits the fact that the ampli- 
tudes for these processes differ in a pure isospin 
amplitude A3/2, defined as the matrix element of 
the isovector part of the effective Hamiltonian be- 
tween a B meson and the ttK isospin eigenstate 
with / = §. In the Standard Model the param- 
eters of this amplitude are determined, up to an 
overall strong phase (f), in the limit of SU(3) fla- 
vor symmetry B. Using the QCD factorization 
theorem proved in Iq], the SU(3)-breaking correc- 
tions can be calculated in a model-independent 
way up to nonfactorizable terms that are power- 
suppressed in A/mi, and vanish in the heavy- 
quark limit. 

A convenient parameterization of the decay 
amplitudes ^+0 = A{B^ -^ 7r+X°) and Ao+ = 
-V2AiB+ -^ 7t"K+) is § 



A 



+0 



P(l-eae*''e"'), 



(1.1) 



Ao+ = P 



1 - e, e^^e'" - £3/2 e'^ 



i4>(„n 



few) 



where P is the dominant penguin amplitude de- 
fined as the sum of all terms in the i?+ -^ Tr~^K° 
amplitude not proportional to e''^ , rj and (j) are 
strong phases, and Ea, £3/2 and Sew are real 
hadronic parameters. The weak phase 7 changes 
sign under a CP transformation, whereas all other 
parameters stay invariant. 

Based on a naive quark-diagram analysis one 
would not expect the 5+ -^ tt^K^ amplitude to 
receive a contribution from b -^ uus tree topolo- 
gies; however, such a contribution can be in- 
duced through final-state rescattering or anni- 
hilation contributions @, 0, |l|, |l|, |||, |2l). 
They are parameterized by £a — 0{\^). In the 
heavy-quark limit this parameter can be calcu- 
lated and is found to be very small, £a ~ —2% 
|22] . In the future, it will be possible to put upper 
and lower bounds on Sa by comparing the CP- 
averaged branching ratios for the decays B^ -^ 
Tr^K° and B± -^ K^K" [|o). Below we assume 
\£a\ < 0.1; however, our results will be almost 
insensitive to this assumption. 



The terms proportional to £3/2 in (1.1) pa- 
rameterize the isospin amplitude A3/2 ■ The weak 
phase e^'^ enters through the tree process b -^ 
uus, whereas the quantity Sew describes the ef- 
fects of electroweak penguins. The parameter 
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£3/2 measures the relative strength of tree and 
QCD penguin contributions. Information about 
it can be derived by using SU(3) flavor symme- 
try to relate the tree contribution to the isospin 
amplitude ^3/2 to the corresponding contribu- 
tion in the decay B^ -^ tt+tt^. Since the final 
state tt+tt'^ has isospin / = 2, the amplitude 
for this process does not receive any contribu- 
tion from QCD penguins. Moreover, electroweak 
penguins in 6 ^ dqq transitions are negligibly 
small. We define a related parameter £3/2 by 
writing £3/2 = £3/2 \/l - 2£q cosr/cos7 + e^^ so 
that the two quantities agree in the limit Ea —^ 0. 
In the SU(3) limit this new parameter can be de- 
termined experimentally form the relation |g[| 



£3/2 



= Ri 



Vus 



Vud 



2B{B^ 



= 7rO) 



-,1/2 



B(B± 



--K0) 



(1.2) 



SU(3)-breaking corrections are described by the 
factor Ri = 1.22 ± 0.05, which can be calculated 
in a model-independent way using the QCD fac- 
torization theorem for nonlcptonic decays p^ . 
The quoted error is an estimate of the theoret- 
ical uncertainty due to corrections oi 0{^^). 
Using preliminary data reported by the CLEG 
Collaboration [E3| to evaluate the ratio of the 
CP-averaged branching ratios in (|l.2|) we obtain 



£3/2 



0.21 ±0.06cxp± 0.01 



th ■ 



(1.3) 



With a better measurement of the branching ra- 
tios the uncertainty in £3/2 will be reduced sig- 
nificantly. 

Finally, the parameter 



<^F,w — R- 



y^Vc. 



v:,Vus 



Xt 



Stt sin 6w 
(0.64 ± 0.09) X ^-^^^ 



1 



31na;t 



Xt 



1 



(1.4) 



Vub/V,b\ ' 

with Xt = {rat/raw)'^, describes the ratio of elec- 
troweak penguin and tree contributions to the 
isospin amplitude A^/2. In the SU(3) limit it is 
calculable in terms of Standard Model parame- 
ters g Q. SU(3)-breaking as well as small elec- 
tromagnetic corrections are accounted for by the 
quantity R2 = 0.92 ± 0.09 0, |2|. The error 
quoted in ( |1.4D includes the uncertainty in the 
top-quark mass. 

Important observables in the study of the 
weak phase 7 are the ratio of the CP-averaged 



branching ratios in the two B 
modes. 



TiK decay 



^*^ 2B(T^T.tfi) - 0-^5 ±0.28, (1.5) 

and a particular combination of the direct CP 
asymmetries, 

Acp(B±^7rOif±) 



A 



Rt, 
= -0.52 ±0.42. 



A I D± 



ICP 



{B- 



(1.6) 



The experimental values of these quantities are 
derived using preliminary CLEO data reported 
in (231 . The theoretical expressions for i?* and A 



obtained using the parameterization in (1.1) are 
R^^ = 1 + 2£3/2 cos (/) {Sew — cos 7) 

+ ^3/2(1 - 2(5ewCOS7 4- (5|w) + 0(£3/2ea) , 

A = 2£3/2sin7sin(/) + 0(£3/2£Q). (1.7) 

Note that the rescattering effects described by 
£a are suppressed by a factor of £3/2 and thus 
reduced to the percent level. Explicit expressions 
for these contributions can be found in |^. 

2. Lower bound on 7 and constraint 
in the (p, fj) plane 

There are several strategies for exploiting the 
above relations. From a measurement of the ra- 
tio R^ alone a bound on cos 7 can be derived, im- 
plying a nontrivial constraint on the Wolfcnstcin 
parameters p and fj defining the apex of the uni- 
tarity triangle M. Only CP-averaged branching 
ratios are needed for this purpose. Varying the 
strong phases 4> and 77 independently we first ob- 
tain an upper bound on the inverse of i?* . Keep- 
ing terms of linear order in Sa yields 

R*^ < (1 + £3/2|(5ew-cos7|) +£^/2sin27 



+ 2£3/2|£a|sin 7. 



(2.1) 



Provided i?* is significantly smaller than 1, this 
bound implies an exclusion region for cos 7 which 
becomes larger the smaller the values of i?* and 
£3/2 are. It is convenient to consider instead of 
R^ the related quantity ]13[| 



Xr,= 



^3/2 



= 0.72 ±0.98cxp± 0.03th. 



(2.2) 
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Because of the theoretical factor i?i entering the 



definition of £3/2 in (1^) this is, strictly speak- 
ing, not an observable. However, the theoretical 
uncertainty in Xjj is so much smaller than the 
present experimental error that it can be ignored. 
The advantage of presenting our results in terms 
of Xj^ rather than i?» is that the leading depen- 
dence on £3/2 cancels out, leading to the simple 
bound \Xii\ < |<5ew - COS7I + 0(83/2, ea)- 




Pi 
X 



25. 50. 75. 100. 125. 150. 175. 
gamma 

Figure 1: Theoretical upper bound on the ratio Xr 
versus I7I for Sa = 0.1 (solid line) and Ea = (dashed 
line). The horizontal line and band show the current 
experimental value with its la variation. 

In Figure n^ we show the upper bound on Xji 
as a function of I7I , obtained by varying the input 
parameters in the intervals 0.15 < £3/2 < 0.27 
and 0.49 < Sew < 0.79 (corresponding to using 
\Vub/Vcb\ = 0.085±0.015in (O)). Note that the 



effect of the rescattering contribution parameter- 
ized by Ea is very small. The gray band shows the 
current value of Xn, which clearly has too large 
an error to provide any useful information on 7.^ 
The situation may change, however, once a more 
precise measurement of Xn will become avail- 
able. For instance, if the current central value 
Xji = 0.72 were confirmed, it would imply the 
bound I7I > 75°, marking a significant improve- 
ment over the indirect limit I7I > 37° inferred 
from the global analysis of the unitarity triangle 
including information from K-K mixing |j]. 

So far, as in previous work, we have used the 
inequality (2.1) to derive a lower bound on |7|. 



However, a large part of the uncertainty in the 
value of ^EW, and thus in the resulting bound 




Figure 2: Theoretical constraints on the Wolfen- 
stein parameters (p, ff) implied by a measurement 
of the ratio Xr in B^ -^ tvK decays (solid lines), 
semileptonic B decays (dashed circles), and Bd,s- 
Bd,s mixing (dashed-dotted line). 



on I7I, comes from the present large error on 
\Vub\- Since this is not a hadronic uncertainty, it 



is appropriate to separate it and turn (2.1) into a 
constraint on the Wolfenstein parameters p and 
fj. To this end, we use that cos 7 
by definition, and few 



p/^/^ 



(0.24±0.03)/A/p2_^r?2 
from (|l.4| ) . The solid lines in Figure || show the 
resulting constraint in the (p, fj) plane obtained 
for the representative values Xn = 0.5, 0.75, 1.0, 
1.25 (from right to left), which for £3/2 = 0.21 
would correspond to i?, = 0.82, 0.75, 0.68, 0.63, 
respectively. Values to the right of these lines 
are excluded. For comparison, the dashed cir- 
cles show the constraint arising from the mea- 
surement of the ratio iKb/Kb] = 0.085 ± 0.015 
in semileptonic B decays, and the dashed-dotted 
line shows the bound implied by the present ex- 
perimental limit on the mass difference Arris in 
the Bs system [||. Values to the left of this 
line are excluded. It is evident from the figure 
that the bound resulting from a measurement 
of the ratio Xji in B^ -^ ttK decays may be 
very nontrivial and, in particular, may eliminate 
the possibility that 7 = 0. The combination of 
this bound with information from semileptonic 
decays and B-B mixing alone would then deter- 
mine the Wolfenstein parameters p and fj within 
narrow ranges,^ and in the context of the CKM 
model would prove the existence of direct CP vi- 
olation in B decays. 



^Unfortunately, the 2cr deviation from 1 indicated by 
a first preliminary CLEO result has not been confirmed 
by the present data. 



^An observation of CP violation, such as the measure- 
ment of e/f in K-K mixing or sin 2/3 in iJ — > J/ip Kg 
decays, is however needed to fix the sign of fj. 
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3. Extraction of 7 

Ultimately, the goal is of course not only to de- 
rive a bound on 7 but to determine this pa- 
rameter directly from the data. This requires 
to fix the strong phase in (|l.7|), which can be 
achieved either through the measurement of a CP 
asymmetry or with the help of theory. A strat- 
egy for an experimental determination of 7 from 
B^ -^ nK decays has been suggested in p. It 
generalizes a method proposed by Gronau, Ros- 
ner and London |g5[ to include the effects of elec- 
troweak penguins. The approach has later been 
refined to account for rescattering contributions 
to the B^ -^ TT^K^ decay amplitudes 0. Before 
discussing this method, we will first illustrate an 
easier strategy for a theory-guided determination 
of 7 based on the QCD factorization theorem for 
nonleptonic decays |^. This method does not 
require any measurement of a CP asymmetry. 

3.1 Theory-guided determination 

In the previous section the theoretical predictions 
for the nonleptonic B — > ttK decay amplitudes 
obtained using the QCD factorization theorem 
were used in a minimal way, i.e., only to calculate 
the size of the SU(3)-breaking effects parameter- 
ized by i?i and R2 ■ The resulting bound on 7 and 
the corresponding constraint in the (p, fj) plane 
are therefore theoretically very clean. However, 
they are only useful if the value of X^ is found 
to be larger than about 0.5 (see Figure P, in 
which case values of I7I below 65° are excluded. 
If it would turn out that Xfi < 0.5, then it is 



possible to satisfy the inequality (2.1) also for 



small values of 7, however, at the price of having 
a very large strong phase, (j) ~ 180°. But this 
possibility can be discarded based on the model- 
independent prediction that B 



4> = 0[asimb),A/mb\ 



(3.1) 



A direct calculation of this phase to leading power 
in A/mb yields 4> ~ -11° 1^- Using the fact 
that (j) is parametrically small, we can exploit a 
measurement of the ratio Xn to obtain a deter- 
mination of I7I - corresponding to an allowed 
region in the (p, fj) plane - rather than just a 
bound. This determination is unique up to a 



sign. Note that for small values of 4> the im- 
pact of the strong phase in the expression for 



i?, in ( |l.7D is a second-order effect. As long as 



\(j)\ ^ ^2Ae3/2/£3/2) the uncertainty in cosi^ 
has a much smaller effect than the uncertainty 
in £3/2- With the present value of £3/2 this is 
the case as long as \(f>\ <C 43°. We believe it is 
a safe assumption to take \(j)\ < 25° (i.e., more 
than twice the value obtained to leading order in 
A/mfc), so that cos<j) > 0.9. 




Figure 3: Allowed regions in the (p, ry) plane for 
fixed values of Xr, obtained by varying all theoret- 
ical parameters inside their respective ranges of un- 
certainty, as specified in the text. The sign of fj is 
not determined. 

Solving the equation for R^ in (IJ) for cos 7, 
and including the corrections of 0{sa), we find 



cos 7 — S 



Xf 



EW 



2^3/2iXj^ 



1 



-^Iw) 



cos 0-1- £3/2(JEW 



_^^^^cos^7smj_^ (3.2) 

cos (j) + e3/2<5EW 

where we have set cos (j) — I in the numerator 
of the 0{ea) term. Using the QCD factoriza- 
tion theorem one finds that £QC0S7y « —0.02 
in the heavy-quark limit p3|, and we assign a 
100% uncertainty to this estimate. In evaluating 



the result (3.2) we scan the parameters in the 
ranges 0.15 < £3/2 < 0.27, 0.55 < Sew < 0.73, 
-25° <(f>< 25°, and -0.04 < £aCosr/sin27 < 0. 
Figure || shows the allowed regions in the (p, fj) 
plane for the representative values Xj^ — 0.25, 
0.75, and 1.25 (from right to left). We stress 
that with this method a useful constraint on the 
Wolfenstein parameters is obtained for any value 
of Xr. 

3.2 Model-independent determination 

It is important that, once more precise data on 
B^ — > ttK decays will become available, it will 
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25. 50. 75. 100. 125. 150. 175. 
gamma 

Figure 4: Contours of constant _R« ("hyperbolas") 
and constant \A\ ("circles") in the (I7I, |(/f>|) plane. 
The sign of the asymmetry A determines the sign 
of the product sin 7 sin (f>. The contours for _R« refer 
to values from 0.6 to 1.0 in steps of 0.1, those for 
the asymmetry correspond to 5%, 15%, and 25%, as 
indicated. 



be possible to test the prediction of a small strong 
phase </) experimentally. To this end, one must 



determine the CP asymmetry A defined in (1.6) 



in addition to the ratio i?*. From (1/7) it follows 
that for fixed values of £3/2 and Sew the quanti- 
ties i?, and A define contours in the (7, 0) plane, 
whose intersections determine the two phases up 
to possible discrete ambiguities |g, 0. Figure ^ 
shows these contours for some representative val- 
ues, assuming £3/2 = 0.21, S^w = 0.64, and Sa = 
0. In practice, including the uncertainties in the 
values of these parameters changes the contour 
lines into contour bands. Typically, the spread 
of the bands induces an error in the determina- 
tion of 7 of about 10° 0.'^ In the most general 
case there are up to eight discrete solutions for 
the two phases, four of which are related to the 
other four by a sign change {j,4>) -^ (— 7,— 0). 
However, for typical values of i?* it turns out that 
often only four solutions exist, two of which are 
related to the other two by a sign change. The 
theoretical prediction that (j) is small implies that 



''A precise determination of this error requires know- 
ledge of the actual values of the observables. Gronau and 
Pirjol Pq] find a larger error for the special case where 
the product j sin 7 sin 01 is very close to 1, which however 
is highly disfavored because of the expected smallness of 
the strong phase </>■ 



solutions should exist where the contours inter- 
sect close to the lower portion in the plot. Other 
solutions with large are strongly disfavored. 



Note that according to (1.7) the sign of the CP 



asymmetry A fixes the relative sign between the 
two phases 7 and 0. If we trust the theoretical 
prediction that cj) is negative |g2[, it follows that 
in most cases there remains only a unique solu- 
tion for 7, i.e., the CP-violating phase 7 can he 
determined without any discrete ambiguity. 

Consider, as an example, the hypothetical 
case where i?* = 0.8 and A — —15%. Figure ^ 
then allows the four solutions where (7,0) ~ 
(±82°,t21°) or (±158°,t78°). The second pair 
of solutions is strongly disfavored because of the 
large values of the strong phase (j). From the first 
pair of solutions, the one with (f) « —21° is clos- 
est to our theoretical expectation that w —11°, 
hence leaving 7 w 82° as the unique solution. 

4. Sensitivity to New Physics 

In the presence of New Physics the theoretical 
description of B^ — > ttK decays becomes more 
complicated. In particular, new CP-violating con- 
tributions to the decay amplitudes may be in- 
duced. A detailed analysis of such effects has 
been presented in p3|. A convenient and com- 
pletely general parameterization of the two am- 



plitudes in (1.1) is obtained by replacing 



P -^ P' . 



ipe 



itpp 



JEW 



ae 



^(pa 



+ ibe 



iipb 



(4.1) 



where p, a, b are real hadronic parameters, and 
4>p, (pa, </*& are strong phases. The terms ip and 
ib change sign under a CP transformation. New 
Physics effects parameterized by P' and p are 
isospin conserving, while those described by a 
and b violate isospin symmetry. Note that the 
parameter P' cancels in all ratios of branching 
ratios and thus does not affect the quantities i?* 
and Xfi as well as any CP asymmetry. Because 



the ratio i?* in (1^) would be 1 in the limit of 
isospin symmetry, it is particularly sensitive to 
isospin-violating New Physics contributions. 

New Physics can affect the bound on 7 de- 
rived from (^j) as well as the extraction of 7 
using the strategies discussed above. We will dis- 
cuss these two possibilities in turn. 
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4.1 EfTects on the bound on 7 



The upper bound on R^^ in (2.1) and the cor- 
responding bound on Xji shown in Figure H are 
model-independent results vahd in the Standard 
Model. Note that the extremal value of R~^ is 
such that \Xfi\ < (1 -I- ^ew) irrespective of 7. A 
value of \Xr\ exceeding this bound would be a 
clear signal for New Physics [0, |lO[ p^ . 

Consider first the case where New Physics 
may induce arbitrary CP-violating contributions 
to the B — > ttK decay amplitudes, while preserv- 
ing isospin symmetry. Then the only change with 
respect to the Standard Model is that the param- 
eter p may no longer be as small as 0{ea)- Vary- 
ing the strong phases and (j)p independently, 
and allowing for an arbitrarily large New Physics 
contribution to p, one can derive the bound U3 



\Xb.\ < Y 1 - 2^EW C0S7 + 5|^ < 1 + Sew ■ 

(4.2) 
Note that the extremal value is the same as in 
the Standard Model, i.e., isospin-conserving New 
Physics effects cannot lead to a value of \Xii\ 
exceeding (1 -|- ^ew)- For intermediate values of 
7 the Standard Model bound on Xji is weakened; 
but even for large values p — 0(1), corresponding 
to a significant New Physics contribution to the 
decay amplitudes, the effect is small. 

If both isospin-violating and isospin-conser- 
ving New Physics contributions are present and 
involve new CP-violating phases, the analysis be- 
comes more complicated. Still, it is possible to 
derive model-independent bounds on Xj^. Allow- 
ing for arbitrary values of p and all strong phases, 
one obtains ll31 



\Xr\ < v/(|a| + |cos7|)2-)-(|5| + |sin7|)2 

< l + ^/a^+b^<—+XR, (4.3) 

^3/2 

where the last inequality is relevant only in cases 
where ^/oF+W ^ 1. The important point to 
note is that with isospin-violating New Physics 
contributions the value of \Xji\ can exceed the 
upper bound in the Standard Model by a poten- 
tially large amount. For instance, if Va^ + b^ 
is twice as large as in the Standard Model, cor- 
responding to a New Physics contribution to the 
decay amplitudes of only 10-15%, then \X]i\ could 



be as large as 2.6 as compared with the maximal 
value 1.8 allowed in the Standard Model. Also, in 
the most general case where b and p are nonzero, 
the maximal value \Xfi\ can take is no longer 
restricted to occur at the endpoints 7 = 0° or 
180°, which are disfavored by the global analysis 
of the unitarity triangle Q. Rather, \X[i\ would 
take its maximal value if | tan7| = \p\ = \b/a\. 



The present experimental value of X^j in (2.2) 
has too large an error to determine whether there 
is any deviation from the Standard Model. If Xn 
turns out to be larger than 1 (i.e., at least one 
third of a standard deviation above its current 
central value) , then an interpretation of this re- 
sult in the Standard Model would require a large 
value I7I > 91° (see Figure P, which would be 
difficult to accommodate in view of the upper 
bound implied by the experimental constraint on 
Bg-Bs mixing, thus providing evidence for New 
Physics. If Xr > 1.3, one could go a step further 
and conclude that the New Physics must neces- 



sarily violate isospin |13 



4.2 Effects on the determination of 7 

A value of the observable i?* violating the bound 



(2.1) would be an exciting hint for New Physics. 
However, even if a future precise measurement 
will give a value that is consistent with the Stan- 
dard Model bound, B^ -^ nK decays provide an 
excellent testing ground for physics beyond the 
Standard Model. This is so because New Physics 
may cause a significant shift in the value of 7 
extracted using the strategies discussed in Sec- 
tion pi leading to inconsistencies when this value 
is compared with other determinations of 7. 

A global fit of the unitarity triangle combin- 
ing information from semileptonic B decays, B- 
B mixing, CP violation in the kaon system, and 
mixing-induced CP violation in B — > J/iJj Ks de- 
cays provides information on 7 which in a few 
years will determine its value within a rather nar- 
row range y]. Such an indirect determination 
could be complemented by direct measurements 
of 7 using, e.g., B -^ DK'^*'> decays (see below), 
or using the triangle relation 7 = 180° — a — /3 
combined with a measurement of a. We will as- 
sume that a discrepancy of more than 25° be- 
tween the "true" 7 — arg(V^*^) and the value JttK 
extracted in B^ — > ttK decays will be observable 
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after a few years of operation at the B factories. 
This sets the benchmark for sensitivity to New 
Physics effects. 




25. 50. 75. 100. 125. 150. 175. 
gamma 

Figure 5: Contours of constant Xr versus 7 and the 
parameter a, assuming 7 > 0. The horizontal band 
shows the value of a in the Standard Model. 

In order to illustrate how big an effect New 
Physics could have on the extracted value of 7 we 
consider the simplest case where there are no new 
CP-violating couplings. Then all New Physics 



contributions in (4.1) are parameterized by the 
single parameter onp = a — Sew ■ A more general 
discussion can be found in |l3|. We also assume 
for simplicity that the strong phase 4> is small, 



as suggested by (3.1). In this case the difference 
between the value ^t^k extracted from B^ — > nK 
decays and the "true" value of 7 is to a good 
approximation given by 

cos 7^i<- ~ cos 7 — flNP . (4.4) 

In Figure ^ we show contours of constant Xf( ver- 
sus 7 and a, assuming without loss of generality 
that 7 > 0. Obviously, even a moderate New 
Physics contribution to the parameter a can in- 
duce a large shift in 7. Note that the present 
central value of Xn « 0.7 is such that values of 
a less than the Standard Model result a « 0.64 
are disfavored, since they would require values 
of 7 exceeding 100°, in conflict with the global 
analysis of the unitarity triangle M . 

4.3 Survey of New Physics models 

In |l3], we have explored how physics beyond 
the Standard Model could affect purely hadronic 
FCNC transitions of the type b -^ sqq focusing, 
in particular, on isospin violation. Unlike in the 
Standard Model, where isospin-violating effects 



in these processes are suppressed by electroweak 
gauge couplings or small CKM matrix elements, 
in many New Physics scenarios these effects are 
not parametrically suppressed relative to isospin- 
conserving FCNC processes. In the language of 
effective weak Hamiltonians this implies that the 
Wilson coefficients of QCD and electroweak pen- 
guin operators are of a similar magnitude. For a 
large class of New Physics models we found that 
the coefficients of the electroweak penguin opera- 
tors are, in fact, due to "trojan" penguins, which 
are neither related to penguin diagrams nor of 
electroweak origin. 

Specifically, we have considered: (a) mod- 
els with tree-level FCNC couplings of the Z bo- 
son, extended gauge models with an extra Z' 
boson, supersymmetric models with broken R- 
parity; (b) supersymmetric models with R-parity 
conservation; (c) two-Higgs-doublet models, and 
models with anomalous gauge-boson couplings. 
Some of these models have also been investigated 
in [|ri|, |l2[. In case (a), the resulting electroweak 
penguin coefficients can be much larger than in 
the Standard Model because they are due to tree- 
level processes. In case (b), these coefficients can 
compete with the ones of the Standard Model 
because they arise from strong-interaction box 
diagrams, which scale relative to the Standard 
Model hke {as/a){ml^,'/■m'g^gY). In models (c), 
on the other hand, isospin-violating New Physics 
effects are not parametrically enhanced and are 
generally smaller than in the Standard Model. 

For each New Physics model we have ex- 
plored which region of parameter space can be 
probed by the B^ -^ nK observables, and how 
big a departure from the Standard Model pre- 
dictions one can expect under realistic circum- 
stances, taking into account all constraints on 
the model parameters implied by other processes. 
Table summarizes our estimates of the maximal 
isospin-violating contributions to the decay am- 
plitudes, as parameterized by |aNp|. They are 
the potentially most important source of New 
Physics effects in B^ -^ ttK decays. For com- 
parison, we recall that in the Standard Model 
a « 0.64. Also shown are the corresponding max- 
imal values of the difference IjttK — j\- As noted 
above, in models with tree-level FCNC couplings 
New Physics effects can be dramatic, whereas in 
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Model 


IanpI \i-^k-i\ 


FCNC Z exchange 

extra Z' boson 

SUSY without R-parity 


2.0 180° 
14* 180° 
14* 180° 


SUSY with R-parity 


0.4 25° 
1.3 180° 


two-Higgs-doublet mod. 
anom. gauge-boson coupl. 


0.15 10° 
0.3 20° 



Table 1: Maximal contributions to onp and shifts 
in 7 in extensions of the Standard Model. En- 
tries marked with "*" are upper bounds derived us- 
ing (4.3). For the case of supersymmetric models 



with R-parity the first (second) row corresponds to 
maximal right-handed (left-handed) strange-bottom 
squark mixing. For the two-Higgs-doublet models 
we take m[j+ > 100 GeV and tan/3 > 1. 

supersymmetric models with R-parity conserva- 
tion isospin-violating loop effects can be compet- 
itive with the Standard Model. In the case of 
supersymmetric models with R-parity violation 
the bound ( [1.3| ) implies interesting limits on cer- 
tain combinations of the trilinear couplings A^ j, 
and A"- J., as discussed in |1^ . 

5. Alternative approaches and recent 
developments 

We will now review recent developments regard- 
ing other approaches to determining 7, focusing 
mainly on proposals that can be pursued at the 
first-generation B factories. 

5.1 Variants of the B^ — > tvK strategy 

The first proposal to constraining 7 using CP- 
averaged B —^ nK branching ratios was made 
by Fleischer and Manncl [pTl, who suggested to 
consider the ratio 



R 



t{B+) B(S0 -^ t:^K^ 



t{B") B(B± ^ 7r±i^0) 



1.11 ±0.33. 



(5.1) 

Neglecting the small rescattering contribution to 
the B^ -^ TT^K^ decay amplitudes as well as 
electroweak penguin contributions yields 

R ~ 1 — 2eT cos <j)T cos 7 -I- £3- 

= sin^7 -I- (et — cos 0t cos 7)^ + sin^^T cos^7 



where St is a real parameter of similar magni- 
tude as £3/2, and (fir is a strong phase. If the 
ratio R was found significantly less than 1, the 
above inequality would imply an exclusion region 
around 7 = 90°. 

Unlike the parameter £3/2 in B — > ttK de- 
cays, the quantity Et is not constrained by SU(3) 
symmetry and cannot be determined experimen- 
tally. The strategy proposed in J2^ is to elimi- 
nate this quantity in deriving a bound on 7. This 
weakens the handle on the weak phase except for 
the particular case where Et « cos0t cos 7. The 
neglect of electroweak penguin and rescattering 
corrections is questionable and has given rise to 
some criticism jl^, !§, |^, |l|. Yet, although 



the bound (5.2) is theoretically not as clean as 



the corresponding bound (2.1) on the ratio i?,, 
a precise measurement of the ratio R can pro- 
vide for an interesting consistency check. Var- 
ious refinements and extensions of the original 
Fleischer-Mannel strategy are discussed in P8[ . 
Some authors have suggested to eliminate 
the small rescattering contribution to the B^ -^ 
TT K^ decay amplitudes, parameterized by £a in 



(1.1), by assuming SU(3) symmetry and exploit- 
ing amplitude relations connecting B — > tiK and 
B — > TTTT decays with other decay modes, such 
as B+ -^ K+K^, B+ -^ Tr+TyW, Bs ^ Kn ^, 
or B+ -> K+rj^'\ B, -> n^rj^'^ ^. Note that 
approaches based on Bg decays have to await 
second-generation B factories at hadron collid- 
ers. Besides relying on the assumption of SU(3) 
flavor symmetry the above proposals suffer from 
theoretical uncertainties related to rj-rj' mixing. 
Given that the rescattering contribution in ques- 
tion is expected to be very small |g^, and that 
this expectation can be tested experimentally us- 
ing B^ -^ K^K^ decays |20 , neglecting Sa or 



> sin 7 , 



(5.2) 



putting an upper bound on it will most likely be 
a better approximation than neglecting the po- 
tentially much larger SU(3)-breaking corrections 
in the above strategies. 

5.2 SU(3) relations between Bd and Bg de- 
cay amplitudes 

Fleischer has recently suggested to use the U-spin 
{d ^^ s) subgroup of flavor SU(3) to derive rela- 
tions between various decays into CP eigenstates. 
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such as fplf 



Bd^J/i^Ks 



B, 



IT TT 



Bs 
Bs 
B, 



DtD: , 

K+K- . (5.3) 



The first two relations provide access to the weak 
phase 7, while the third one is sensitive to both 
(3 and 7. Although this strategy involves B^ de- 
cays and thus cannot be pursued at the first- 
generation B factories, we discuss it because of 
its general nature. 

Consider the example of Bd,Bs — > J/tpKs 
decays, which are governed by an interference of 
tree and penguin topologies. The sensitivity to 
7 arises from the presence of top- and up-quark 
penguins. A general parameterization of the de- 
cay amplitudes Ad,s = A{Bd,s -^ J/"^ Ks) is 



Ad = A'{l + tan^ece'e"*' e 



' pi4>' pil\ 



As ^ Ail 



ai<P c«7 



), 



where 9c is the Cabibbo angle, </>'■'' are strong 
phases, and the penguin contributions are pro- 
portional to parameters e^^'^ ^ 0.2. Exact U- 
spin symmetry would imply A' = A, e' = e, 
and (j)' = 0. In that limit 7 could be deter- 
mined (with discrete ambiguities) by measuring 
the direct and mixing-induced CP asymmetries 
A^§.iB, ^ J/^Ks) and A^^iB, ^ J/^Ks) 
as well as the ratio of the CP-averaged Bd,Bs -^ 
J/ip Ks decay rates. These observables would 
suffice to fix e, (j) and 7. Assuming A' — A, this 
parameter cancels out. 

This approach is interesting in that it is gen- 
eral and can be applied to several different decay 
modes |pl|. Assuming the theoretical uncertain- 
ties related to U-spin breaking can be controlled, 
it will provide several independent determina- 
tions of (3 and 7. However, a question that needs 
to be addressed in future work is how important 
the SU(3)-breaking corrections leading io A' ^ A 
are, given that e ~ 0.2 is expected to be a small 
parameter. 

5.3 Dalitz-plot analysis in B^ -^ ■k^tv~^tv~ 
decays 

There have been several proposals for obtaining 
information on the weak phases a and 7 from 



an analysis of the Dalitz plot in B ^- Stt decays. 
In fact, the approach of Quinn and Snyder |32|1 
based on the decays _B — > /ott — > 37r is considered 
to offer one of the most promising ways to mea- 
sure a [ p3[. The strategy of Bediaga et al. Q 
(see also [p5[) for learning 7 is to fit the measured 
Dalitz distribution d F/dm^ dm^ in the decays 



B^ — > TT+TT+TT , where mf 
an ansatz of the form 



(P^+ +P7r-) , to 



J2a^e"''F^i^^lml) 



(5.5) 



Fi{ra\,m\) are appropriate kinematic functions 
for resonance or continuum contributions. From 
the fit one extracts a set {ai,9i\ of real ampli- 
tudes and complex phases. Performing a similar 
fit to the CP-conjugated decays B~ — > 7r~7r~7r+ 
gives parameters {ai,Oi}. The complex phases 
are sums of strong and weak phases, and the 
latter ones are determined from the differences 



(5.4) 0, 



i — Oi). The weak phase 7 enters through 



the interference of the Xc"" resonance state with 
other resonance channels (e.g., p'^n^, /ott^, etc.) 
and nonresonant contributions. The associated 
CKM factors are V*ijVcd and V*i^Vud, respectively. 

A theoretical problem inherent in this ap- 
proach is the "penguin pollution", i.e., the fact 
that b — > dqq penguin transitions contaminate 
the analysis. If the penguin/tree ratio is assumed 
to be at most 20%, then the resulting error in the 
extraction of 7 is bound to be less than 11° p4| . 
Unfortunately, the recent data on _B — > ttK and 
B ^ TTTT decays reported by the CLEO Collab- 
oration |2^ suggest that the penguin/tree ratio 
may be significantly larger than 20%. 

The feasibility of this method profits from 
the fact that no tagging is required (only charged 
i3-meson decays are considered), the final state 
consists of three charged pions (no tt'^ reconstruc- 
tion is needed), and a Dalitz plot analysis typ- 
ically does not require very large data samples. 
The authors of |3J] estimate that with only 1000 
events one could obtain a resolution of A7 « 20°. 
Potential problems of the approach are that the 
size of the interference term depends on the yet 
unknown B^ —> Xc^^^ ~* tt 7r''"7r~ branching ra- 
tio, and that contamination from nonresonant 
channels may, in the end, require larger data 
samples. 
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5.4 Extracting 7 in iJ — *• DK'^*^ decays 

B -^ DK^*^ decays were originally considered to 
be the "classical" way for determining 7. Later, 
it was realized that this is a very challenging 
route, which poses high demands to experiment 
and theory. We discuss this strategy because it 
provides an extraction of 7 from tree processes 
alone, which is unlikely to be affected much by 
New Physics. 

The original idea of Gronau and Wyler pq| 
(see also ^, |8|) was to use the interference of 
the amplitudes for the decays i?+ -^ D^K^ and 
i?+ -^ D'^K^ occurring if the charm meson is 
detected as a CP eigenstate D? = 75 (-0° + D°)- 

The first decay is due to the quark transition 6. Conclusions 
h -^ cus proportional to V*ijVus, whereas the 
second one is due to the decay b -^ ucs pro- 
portional to V*fycs- The relative phase of these 
two combinations of CKM elements is 7. Ideally, 
one would measure the six rates for the decays 
B+ -^ D^K+, B+ -^ D°K+, B+ ^ D'iK+, 
and their CP conjugates. Then, using isospin tri- 
angles, 7 could be determined in a theoretically 
clean way. 

This strategy is hindered by the fact that 
is it not possible experimentally to determine 
the rate of the doubly Cabibbo-suppressed decay 
B+ -^ D"K+ followed by D° ^ K-Tr+, because 
its combined branching ratio is similar to that of 
the transition i?+ -^ D^K^ followed by the dou- 
bly Cabibbo-suppressed decay l)" — > K^tt^ p^ . 
Several approaches have been suggested to cir- 
cumvent this problem p9, Bfl; however, they are 
challenging from an experimental point of view 
because precise measurements of very small de- 
cay rates are required. 

Recently, some authors have suggested to use 
isospin relations combined with certain dynami- 
cal assumptions (the neglect of annihilation con- 
tributions relative to color-suppressed tree ampli- 
tudes) to ehminate the "difficult" B+ -> D°K+ 
and B^ -^ D^K^ rates in favor of six other 
B,B ^ DK rates ||l|, ||. Unfortunately it is 
difficult to gauge the accuracy of the assumptions 
made in these proposals. 

Considering the various options that have 
been discussed it appears that measuring 7 at the 
first-generation B factories using B — > DK^*^ de- 



cays is very challenging |l:3| , and more demand- 
ing than a determination based on B^ -^ ttK 
decays. On the other hand, to have two indepen- 
dent determinations of 7 from these two classes of 
decays would be extremely important. Whereas 
7 measured in i? ^ _D iiT^*' decays is likely to 
be the "true" phase of the CKM matrix, the an- 
gle JttK determined in B^ -^ ttK decays probes 
loop processes and may easily be affected by New 
Physics. As discussed in Section and summa- 
rized in Table H, comparing the two determina- 
tions would provide a very sensitive probe for 
physics beyond the Standard Model. 



Among the strategies for determining the weak 
phase 7 = a,Tg{V*f^) of the quark mixing ma- 
trix, approaches based on rate measurements in 
B — > ttK decays play an important role and 
have received a lot of attention recently. The 
corresponding branching ratios are "large" and 
the final states are "easy" to detect experimen- 
tally. Using isospin, Fierz and flavor symmetries 
together with the fact that nonleptonic B de- 
cays into two light mesons admit a heavy-quark 
expansion, a largely model-independent descrip- 
tion of certain observables concerning the char- 
ged modes B^ -^ ttK can be obtained. Various 
proposals exist for extracting information on 7 
and on the Wolfenstein parameters p and fj us- 
ing these decays. In the future, this will allow 
us to determine 7 with a theoretical uncertainty 
of about 10°. When combined with an alterna- 
tive measurement of 7 using other decays such 
as S — > DK'^*\ this will provide for a sensitive 
probe of physics beyond the Standard Model. 

Acknowledgments 

It is a pleasure to thank the SLAC Theory Group 
for the hospitality extended to me during my visit 
earlier this year. I am grateful to M. Beneke, 
G. Buchalla, Y. Grossman, A. Kagan, J. Ros- 
ner and C. Sachrajda for collaboration on parts 
of the work reported here. This research was 
supported by the Department of Energy under 
contract DE-AC03-76SF00515. 



11 



Heavy Flavours 8, Southampton, UK, 1999 



Matthias Neubert 



References 



A documentation of this measurement can be 



[3] 
[4] 



[7 



Ji 



M 



[lo; 



m 



fl2- 



^ 



4i4 



ft^ 



[16- 



\a 



m 



[19 



w 



found at tittp://www-cdf .fnal.gov/physics/ 



new/bottom/cdf4855. 

A. Alavi-Harati et al. (KTeV Collaboration), 



Phys. Rev. Lett. 83 (1999) 22 



A documentation of this measurement can be 



[22] M. Beneke, G. Buchalla, M. Neubert and C.T. 
Sachrajda, in preparation. 

[23] R. Poling, Rapporteur's talk presented at the 
19th International Symposium on Lepton and 
Photon Interactions at High Energies, Stanford, 
California, 9-14 August 1999; 
CLEG Collaboration, Conference contributions 
CLEG CONF 99-14 and CLEG CONF 99-16. 



found at tittp://www.cern.ch/NA4£ 



For a review, see: Tfie BaBar Pfiysics Book, 



[24' 
[25' 



P.F. Harison and H.R. Quinn eds., SLAG Re- 
port No. SLAC-R-504 (1998). 



M. Neubert and J.R. Rosner, 



(1998) 403 



Pfiys. Lett. B 441 



[26 



[27' 



M. Neubert and J.R. Rosner, Pfiys. Rev. Lett. 



81 (1998) 507e 



M. Neubert, 



014 



J. High. Energy Ptiys. 9902 (1999) 



M. Beneke, G. Buchalla, M. Neubert and C.T. 



[28 



[29 



Sachrajda, Pliys. Rev. Lett. 83 (1999) 1914. 



M. Neubert and B. Stech. in: Heavy Flavours 



[30 



f2nd Edition). A.J. Buras and M. Lindner eds. 



(World Scientific, Singapore, 1998) pp. 294. 
R. Fleischer and J. Matias, Preprint 



iep-ph/9906274 



[31 



[32 



D GlinnHhirry, B Dutta and A KunHii Phya 



Lett. B 456 (1999) 185 



X.-G. He, C.-L. Hsueh and J.-Q. Shi, Preprint 
tiep-ph/9905296 . 



[33 

[34 



^ — Grossman, — M — Neubert and — A-L — Kagan, 



Preprint hep-ph/9909297 



[35 



A . S . Digho , M , Gronau and J . L . Rosner , Phys . 



Rov. Lett. 70 (1007) d ggg 



N.G. Doahpando and X. G. Ho, Phyo. Rev. Lett. 



[36 



74 (1005) 2 g [Erratum: %h%d. 74 (1005) 4000]. 



B. Blok, M. Gronau and J.L. Rosner, Phys. Rev. 



Lett. 78 (1997) 3999 



M. Gronau and J.L. Rosner, 
(1998) 1130051 



Phys. Rev. D 58 



[37" 
[38 

[39 



A, J, Buras, R, Fleisch er and T, Manuel, \^ud. 



Phys. B 533 (1998) 5 



[40 



J.M. Gerard and J. Wo^rore, Eur. Phye. J. C 7 



fl999) 1 



[41 



D. Delepine et al.. \Phvs. Lett. B 429 (1998) 104 



M. Neubert, 



Phys. Lett. B 424 (1998) 152 



[42' 



A.F. Falk et al.. 



Phys. Rev. D 57 (1998) --489e 



D. Atw ood and A. Soni, \Phys. Rev. D 58 (1998) 



036005. 



[43; 



R. Fleischer, Phys. Lett. B 365 (1996) 399 



M. Gronau, J.L. Rosner and D. London, Phys 



Rev. Lett. 73 (1994) 21 



M. Gronau and D. Pirjol, Preprint 
^iep-ph/9902482 



R. Fleischer and T. Manuel, 



(1998) 2752 



Phys. Rev. D 57 



R. Fleischer, Eur. Phys. J. C 6 (1999) 451 



A.F. Buras and R. Fleischer, Preprint 
^iep-ph/981026C . 



M. Gronau and D. Pirjol, Phys. Lett. B 449 



(1999) 321 



K. Agashe and N.G. Deshpande, Phys. Lett. B 
454 (1999)"359| ; \Phys. Lett. B 451 (1999) 215| . 



R. Fleischer, 



Eur. Phys. J. C 10 (1999) 29£ 



Phys. Lett. B 459 (1999) 30e 



A.E. Snyder and H.R. Quinn, Phys. Rev. D 48 



(1993) 2139 



I. Dunietz, these Proceedings. 

I. Bediaga, R.E. Blanco, C 
Mendez-Galain 



Gobel and R. 



4067 



Phys. Rev. Lett. 81 (1998) 



B. Bajc et al., [Phys. Lett. B 447 (1999) 313|; 
S. Fajfer, R.J. Oakes and T.N. Pham, \Phys 



Rev. D 60 (1999) 054029|. 



M. Gronau and D. Wyler, Phys. Lett. B 265 



(1991) 174 



I. Dunietz, 



Phys. Lett. B 270 (1991) 75 



D. Atwood, G. Eilan, M. Gronau and A. Soni, 



Phys. Lett. B 341 (1995) 371 



D. Atwood, I. Dunietz and A. Soni, Phys. Rev. 



Lett. 78 (1997) 3257 



N. Sinha and R. Sinha, 



(1998) 3706 



Phys. Rev. Lett. 80 



111302 



J.H. Jang and P. Ko, Phys. Rev. D 58 (1998) 



M. Gronau, [Phys. Rev. D 58 (1998) 03730l|; 
M. Gronau and J.L. Rosner, Phys. Lett. B 439 



(1998) 171 



A. Soffer, \Phys. Rev. D 60 (1999) 054034 



12 



